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Abstract 

We perform a microscopic evaluation of nuclear generalized parton distributions (GPDs) for spin-0 nuclei in the 
framework of the Walecka model. We demonstrate that the meson (non-nucleon) degrees of freedom dramatically 
influence nuclear GPDs, which is revealed in the non-trivial and unexpected ^-dependence of Deeply Virtual 
Compton Scattering (DVCS) observables. In particular, we find that the first moment of the nuclear D-term, 
(1a(0) oc A 2 ' 26 , which confirms the earlier prediction of M. Polyakov. We find that in the HERMES kinematics, 
contrary to the free proton case, the nuclear meson degrees of freedom in large nuclei enhance the nuclear DVCS 
amplitude which becomes comparable to the Bethe-Heitler amplitude, and, thus, give the non-trivial A-dependence 
to the DVCS asymmetries: as a function of the atomic number the beam-charge asymmetry increases whereas 
the beam-spin asymmetry decreases slowly. 
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I. INTRODUCTION 



During the last decade generalized parton distributions (GPDs) became a standard tool for the 
parameterization of the nonperturbative hadronic structure in hard processes. Although inaccessible 
directly, GPDs enable us to probe the 3-dimensional structure of the target [ll and to study parton 
correlations in the target. GPDs are intensively studied both theoretically [a I5T 4 IEL Ifil IsL M floL ITU 
Il2l llH . 0, 0| and experimentally 0, 0, 0, 0] . One of the key processes used for measurements of 
GPDs is Deeply Virtual Compton Scattering (DVCS). 

While most of theoretical and experimental analyses of DVCS involve the nucleon target, there 
also exist ambitious projects to study DVCS on nuclei. Such experiments may provide us with valuable 
information on the nuclear forces [20| as well as on the change of nucleon properties in the nuclear 
medium. A first DVCS experiment on neon nuclei was performed at HERMES (DESYJhjj. Another 
measurement of nuclear DVCS is planned at the future Electron Ion Collider (EIC) |15l |21| . 

By analogy with inclusive deep inelastic scattering (DIS) on nuclei one can expect that DVCS on 
nuclei will be sensitive to many nuclear phenomena such as shadowing, antishadowing, EMC-effect and 
Fermi motion. In addition, in the DVCS amplitude there maybe present other new nuclear effects which 
are absent in the imaginary part of the forward virtual photon-nucleon scattering amplitude probed in 
inclusive DIS on nuclei. 



DVCS on different nuclei was a subject of investigation in |22l.23, 24, 2 5j|2jJ 123,123 ■ In the last two 
works nuclear GPDs are expressed in terms of convolution (sum) with GPDs of separate nucleons. This 
assumption is rather natural and is based on the well-known fact that nucleons in the nuclei are weakly 
bound and, thus, in hard processes one can consider nucleus as a collective of quasifree nucleons. However, 
as we will show, this assumption is not universal, i.e. it does not work for some of the observables. 

A rather interesting property of nuclear GPDs was predicted in [2(J • Contrary to the expectations 
based on the quasifree nucleon model, the first moment of the D-term, which is intrinsically related 
to GPDs, has nonlinear A-dependence, gU(0) oc A 7 / 3 . Inspired by the result of [2£j], we performed a 
microscopic study of the nuclear GPDs in the framework of the well-established nuclear structure model 
developed by Walecka and collaborators 29, 30, 31]. We confirmed the result of Ref. |2J]j and found rapid 
^4-dependence, d^(0) oc A 2 26 , largely due to mesonic degrees of freedom. We also found that mesons 
significantly enhance the DVCS amplitude compared to the nucleonic contribution. For large nuclei in 
the current HERMES kinematics [l!j, the DVCS amplitude squared increases as |^4wcs| 2 x ^4 4 ' 29 - We 
found that the nuclear asymmetries are very sensitive to the meson internal structure and thus in the 
study of the nuclear DVCS, one should pay particular attention to the meson distributions in the nuclei 
as well as to the quark distributions inside the mesons. We predict that, as a function of the atomic 
number, the beam-charge asymmetry grows as oc A 5 , whereas the beam-spin asymmetry is a slowly 
decreasing function of the atomic number A: Am oc ^4~ 03 . In the absence of mesonic (non-nucleonic) 
degrees of freedom, the asymmetries are virtually independent of A. Also, both asymmetries should 
have a maximum when the DVCS and BH amplitudes have comparable magnitudes. The position of the 
maximum is very sensitive to the nuclear constituents' model. In the forward limit, our GPDs reproduce 
the earlier results for nuclear light-cone distributions [3^, |23, . 

The paper is organized as follows. In the Sect. [H] we demonstrate that the first moment of the 
D-term evaluated in the Walecka model depends on the atomic number A as gU(0) oc A 2 - 26 . In Sect. IIIII 
we evaluate the nucleonic and mesonic off-forward light-cone distributions in nuclei. We discuss the 
influence of mesonic degrees of freedom on physical observables such as the ratio of nuclear-to-nucleon 
GPDs, beam-charge and beam-spin DVCS asymmetries in Sect. II VI In Sect. we summarize our results 
and draw conclusions. 
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II. d A (0) AND ENERGY-MOMENTUM TENSOR 

In this section we perform a microscopic calculation of the ^-dependence of the first moment of the 
D-term, g?a(0), which is intrinsically related to nuclear GPDs. In our analysis we use the connection of 
(1a(0) to the form factors of the energy- momentum tensor introduced in 20] . As a framework we use the 
field-theoretical Walecka model. In this section nuclear constituents (nucleons and mesons) are treated as 
elementary pointlike objects. The influence of their internal structure is considered in the next sections. 

The Lagrangian of the model in the simplest form is [30l l3l| 

C = ${%d — M — g v V + g.flr/, + ^(d^d^ - m 2 </> 2 ) - I V^W" + ^V^V , (1) 

where V^„ — d^V v —d v V^] ip corresponds to the field of nucleons; the massive vector field V and scalar field 
are effective fields, which represent the empirically observed dominant vector and scalar components 
of the NN interaction. The relativistic MFT models based on Lagrangians of type QJ succeeded in 
description of such important characteristics as nuclear densities, the level structure of the nuclear shell 
model, the spin dependence of nucleon-nucleus scattering etc. The simplest version of the model used in 
the work consists of baryons and isoscalar scalar and vector mesons. A virtue of such a simple model is that 
it has only two independent parameters (coupling constants) which can be fixed from the infinite nuclear 
matter properties. After this, all the model predictions for finite nuclei become unambiguous. Extensions 
of the model with more independent degrees of freedom (and consequently more free parameters) can 
give better quantitative description of nucleus 53]. Note that we neglect the pseudoscalar pion degrees of 
freedom since the effects of pions in the ground-state of spin-zero nuclei essentially average to zero due to 
the spin-dependence of the pion-nucleon coupling constant. In addition to the pions, in order to achieve 
a truly quantitative description of nuclear properties, one should ultimately include the rho- meson and 
electromagnetic fields. However, those will lead to small corrections, which we neglect in our exploratory 
study. The values of the coupling constants are fixed from phenomenological parameters of the nuclear 
matter |3(ij ]. Numerically the constants are large. Therefore, the perturbative methods cannot be used. 
Instead, observing that for sufficiently large nuclei, the nuclear density becomes large, one can use the 
mean-field approximation, when the quantum meson fields are replaced by their classical ground-state 
expectation values. The model with Lagrangian (JJJ is not renormalizable and should be understood as 
an effective one. 

Solving classical equations of motion in the nuclear rest frame |35| and making Lorentz boosts to 
the infinite momentum frame, we can obtain the light-cone wave functions necessary for evaluation of 
light-cone matrix elements. 

As it has been shown in 0, one can parameterize the matrix element of the energy- momentum 
tensor in terms of two form factors ild(i) and d A {t)' 

(P'\f^(0)\P) = M A {t)P^P u + ^zd,A {t) (A^ A„ - g^A 2 ) , (2) 

D 

where P = (P + P') /2, A = P' -P. 

Acting with the traceless operator 

\dA* dAi 3 dA k dA k ) Ia=0 < p=const 

on Ecj. J3J) , one can express cIa (0) in terms of the matrix element of the energy-momentum tensor T^ v 
as 20] 

dA(0) = / d 3 r (r irj - S fr^ 2-(r; 0), (3) 
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where m A is the mass of the nucleus; Ty(r;A) is a shorthand notation for the Fourier of the matrix 
element 

1 f r7 3 A 

2E J ^3^> + A/2|T 4 ,(0)|F- A/2) (4) 

between the nuclear ground-states. One can split the energy-momentum tensor, which can be derived 
from the Lagrangian 0), into three parts 

- (-V^Vx„ + ^VxpV^ + rrfyVuVu - g^rn 2 v V p V p ) + (d^d^ - ^{d P 4>d P 4> - m 2 2 )) 

+ (\ + Jud^tp) - gvVfj.^'ipJ , (5) 

which will be referred to as TY^ V ,T^ V and respectively. In complete analogy with J2J one can define 
the form factors Mi/ A (t), di/ A (t) for each of the three contributions T* to the total 7)^. The form factors 
^i/jiWi^/jiW are connected with M^i), dji(i) as 

Afx(t) = J! M i/A{t) = M N/A {t) + M v/A {t) + M HA (t), 

i 

d A{t) = ^2 d l/A (t) = d N/A (t) + d v/A (t) + d^/ A (t). (6) 

i 

In Sect. 11111 wc shall show that the form factors Mi/ A (t), di/ A (t) are closely related to the moments of 
GPDs of nucleons and mesons. Assuming spherical symmetry of the considered nuclei, the final answer 
for d i/A (0) is 

d N /A(0) =-^jd 3 r (r irj - S fr 2 ^j £ $„(r) 7 ^$„(r) - J?- f d 3 r p B (r)V (r), 
«Wa(0) = "T / >rrV'(r) 2 , 

^(0) = ^/rf 3 rrV '(rf, (7) 

where the sum over n in the first equation goes over the occupied levels; 4>„(r), 4>(r), Vo(r) are the 
(classical) self-consistent solutions of the Hartree-Fock equations; ps(r) = ^ n $^(r)<I> n (r) is the baryon 
density. 

The numerical evaluation of d A (0) for different nuclei according to Eq. (JJJ gives the results pre- 
sented in Table H] and Fig. ^ The most interesting feature of the obtained result is that d A (0) receives 
the dominant contribution from the mesons. Notice that <i™ es (0) = gL/^(0) + dy/ A (0) and d A (0) for all 
nuclei lie on the straight lines when plotted in the logarithmic coordinates. The only exception is d A (0) 
for 12 C , where we observe compensation of the nucleon and meson contributions. The least-square fit 
to the values presented in Table [I] gives 

d A (0)^ -0.308 A 2 26 . (8) 

Note that the obtained parameters contain uncertainties due to the numerical nature of our analysis and 
the fitting method. One can see that this result agrees with the estimate based on the liquid drop model 



d A (0) cx -0.2 A 7/3 (l 



A 2/3 



(9) 
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Nucleus 


dT s (0) 


d A (0) 


12 C 


-100.9 


-7.77 


16 


-189.0 


-143.8 


40 Ca 


-1622.5 


-1525 


90 Zr 


-11270.2 


-8258 


208 pb 


-64187.3 


-49195 



TABLE I: Values of gU(0) for different nuclei. Contribution of mesons d™ es (0) = <2<£/a(0) + d v /A{0) and total 
result g(a(0). 



where the last 0(l/j4 2 / 3 )-term takes into account finite width of the nuclear border. Thus the simple 
power dependence ^^(0) ~ A n exists only for relatively large- A nuclei. In our analysis we cannot separate 
contributions of different parts of the nucleus. This explains why the point 12 C on the Fig.QJ does not 
lie on the straight line. 

The A-dependence of the D-term was also examined in [3^ , where it was found that 

d A (0) oc A(l + 0(lnA)). (10) 

As it is discussed in Sect. lIIll the result of [3|J, which is inconsistent with |2(j and the result of our work, 
is due to the nonrelativistic approximation used by the authors. 

It was emphasized in |2(j that the nuclear surface responsible for stability of the liquid drop gives 
the dominant negative contribution to cZyi(0) with rapid A-dependence. In the Walecka model mesons 
compensate the positive contribution of the nucleons to the total pressure and, thus, provide stability of 
the nucleus. A similar phenomenon happens in the calculation of the nucleon £)-term in the chiral quark 
soliton model, where the dominant negative contribution does also come not from the valent quarks but 
from the Dirac sea j23| ■ Thus in all three models we obtain a negative value of (1a (0) which comes from 
"complementary" degrees of freedom. 

To measure D-term, one should study hard processes sensitive to the real part of the DVCS 
amplitude, such as the beam-charge asymmetry 13a o r the DVCS total cross section in properly chosen 
kinematics, where the BH amplitude is suppressed [33. We address this issue in Sect. II VI and show that 
the A-dependence of the GPDs is observable in DVCS asymmetries. 



III. EVALUATION OF NUCLEAR GPDS 



In dealing with hard processes on nuclear targets one often assumes that the hard scattering occurs 
on quasifree nuclear constituents. For the weakly-bound system it also seems quite natural to assume that 
one can ignore possible " distortion" of the constituents in the nuclear medium and neglect offshellness of 
the nucleons and mesons. While this assumption is commonly believed to be justified for the nucleons, 
we cannot estimate its accuracy for the meson degrees of freedom. However, since the meson GPDs are 
unknown even for the onshell case, it should be a sufficient first approximation to use the parametrization 
which only satisfies general GPDs properties. A straightforward application of these ideas to the DVCS 
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FIG. 1: <1a(0) as a function of A. We can see that (1a(0) is mostly due to the meson contribution. 

process results in the convolution formula (we consider only the quark nuclear GPDs) 

H q/A (x,^t) = jH l/A (yA,t)H q/l i,tj , (11) 

where i labels nuclear constituents (ip, V and 4> in case of the Walecka model); Hi/ a describe the dis- 
tribution of the constituents in the nucleus; H q u are the GPDs of the free constituents. (For a more 
detailed discussion, see e.g. |4jJ, ED and references therein. A generalization to the off-forward case and 
non-nucleonic constituents is straightforward.) Equation complies with the intuitive picture of the 
nuclear hard scattering as a two-step process depicted in Fig. the hard scattering amplitude on the 
nuclear target equals the hard amplitude on the free constituent convoluted with the off-forward distri- 
bution of the constituent in the nuclear target. Note that the convolution approximation of Eq. Hll|) 
neglects the simultaneous coherent interaction of the virtual photon with several nuclear constituents. 
Therefore, strictly speaking, Eq. is applicable only for x > 0.1. 

Assuming that the polynomiality conditions 




x n dxH i/A ( X> U)= E 




FIG. 2: Intuitive "Handbag" diagram of nuclear DVCS. 



pi N 

/ X n dx H q/i ( X> Z,t)= ]T C n(t)^ k , 
J — 1 i — n o 



(12) 



fc=0,2,... 



fulfil separately for the constituents' GPDs in the nuclear medium and off-forward nuclear distributions, 
from Eq. Qllf) one can obtain that the polynomiality condition for the total GPD reads 

r i M 

/ x n dxH q/A (x,Z,t)= Yl ^WC fc , (13) 

fe=0,2,... 

where 

[n] 

D k n = Y Cl n B n-l- (14) 
1=0 

For the case of the zero-order moment, summation with the quark charges e q of Eq. (|14|l gives the nuclear 
electromagnetic form factor 

F em (t) = Y e i F N/A{t)F q/N {t) = F N/A (t)F N (t), (15) 
q 

where F/v(i) is the nucleon electromagnetic form factor. Since mesons are not charged, they do not 
contribute to the nuclear electromagnetic form factor. Mathematically it follows from the antisymmetry 
of the mesonic off-forward distributions with respect to the variable x, as it will be shown below. 
From Eq. iQ , one obtains the first moment 



xdxH q/A {x,^t) = M q/A {t) + -Z 2 d q/A (t), (16) 
l o 



where 

M q/A (t)=Y,M q /i(t)M i/A (t), 
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<W*) = +d q /i{t) J ^H i/A (y,£,t)^ , (17) 

and we introduced the conventional notation Mi(t),di(t) for the form factors Bi(t), C\(t). These form 
factors already appeared in the Lorenz-covariant expansion of the energy-momentum tensor in the previ- 
ous section. Note that in our numerical analysis, we use a simple parameterization of the meson GPDs, 
see Eqs. <|21|) an d which corresponds to the vanishing d q j^ and d q / v . 

Contrary to naive expectations, as it follows from Eq. 117fl . the total dA(t) is not reducible to a 
mere sum of the free constituents' cij/^ (i)-terms but consists of two parts. The first term in (|17fl is a 
"collective" effect: in the previous section it was shown that di/A(0) oc A 2 26 . To evaluate the last term 
in (|17fl one can use the approximation H N /A{y,Qify ~ AS (y — 4) which corresponds to an ensemble 
of quasifree nucleons and is justified for a weakly bound nucleus. Then the last term in i|17|) is oc A 2 . 
A direct evaluation with the nucleon off-forward distribution obtained in the Walecka model gives the 
same A-dependence. Thus the predicted A-dependence for the total gU(0) is due to the first term in 
Eq. (|17J) . which dominates in the large- A limit and defines asymptotics, whereas the second term is a 
O ( J 4 _0 - 26 )-correction. Notice also that the essential contribution to the coefficient in front of A 2 26 comes 
from the mesonic degrees of freedom. 

Now we evaluate the nuclear part of the DVCS amplitude Hi/^ix, £, t). In complete analogy with 
Ref. one can define off- forward distributions of nucleons and mesons in the nucleus: 

"v^t) = j-L. (f) v: (-f) |p> 



m 2 , f dz 
' AxP+ 



The definitions of the nucleon and meson GPDs in Eq. I|23|) are chosen such that there is one-to-one 
correspondence to the H — h component of the energy momentum tensor of the model, see Eq. [S] This 
automatically guarantees conservation of the momentum sum rule 



E 



i 

dxxH i/A (x,0,0) = 1. (19) 



Notice that the large-A nucleus in the Mean Field Theory approach in several respects resembles 
the nucleon in the large- N c chiral quark soliton model. A large number of particles in both models enables 
us to represent the interaction of a single particle with the ensemble of other particles as the interaction 
with the central potential. Forward distributions of nucleons (quarks) in both cases are strongly peaked 
at 1/A (1/N C ). So we expect that the off-forward distributions should have similar qualitative features. 

Using standard steps [ITll43| . see also the Appendix for more details, we can obtain from Eq. I|18|) 
H N/A {x^t) = 2£ / |- e — / d*X^n (l-X) P ef-> W - d \ + Z n (-| -X) .(20) 

To simplify evaluations, numerically it seems reasonable to use the approximation V^(x) ~ S^qV in the 
path exponent since actually we have "cutoff' due to the wave functions <&„(?'). We define the average 
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value V as 

V~±Jd 3 XV + (X)p B (X), (21) 

which ensures that the approximate and exact formulas give the same results for the first moment of 
Hn/a (conservation of the nuclear momentum sum rule). This approximation implies that the Wilson 
link just results in a shift of the whole parton distribution by the value Sx oc V . The evaluation with 
the exact expression <|2(Jll shows that Wilson links results not only in the shift but also in the broadening 
of the nucleon distribution. However such "broadening" is negligible, especially for large nuclei. The 
expressions for off-forward distributions look more compact in the momentum representation: 

H N / A {x,£„t) = 

I dJ WT 6 { X + Sx -^)^ # " ( ( * " ^ + ¥) ^ { {x + ^ - ^ 

I ^ (■ - &) ** (<* - ^ - ^ 

5 [ X fTT ) 1-^7 V+ \ {X - £),p± + — )V+ \{x + 0,P±- - -7T -( 22 ) 



(2tt) 3 V P + ) 4.tP+ + V 2 / V ' 2 

We can see that similarly to the gluonic GPD, the vector off- forward distribution is singular at 
x = 0. Using explicit expressions for the off-forward distributions, one can check that they fulfil required 
sum rules. For instance, the momentum sum rule reads: 



x dx H n /a(x, 0, 0) + / x dx H^j^x, 0, 0) + / x dx H v /^(x, 0, 0) = 

= -L J d?X (-V +X (X)V^(X) + m 2 v V + (X)V+(X) + d+<p{x)d+<p{x) 

fd 3 XT++ 

+^{X) 1+ id + ^(X) - g v i,(X) 7+ V + (X)Tp(X)) = = 1 , (23) 

where T ++ is the H — h- component of the energy-momentum tensor and the integral is performed over 
the hypersurface X + = const [3(j. One can easily check that the coefficients in front of £ 2 in the first 
moments J x dx Hi/j\(x, £, t) of the off-forward distributions coincide with the results obtained in the 
previous section from the energy- momentum tensor form factors, see Eq. ( |7| ). 

Plots of the off-forward distributions Hit a as functions of the variable Ax a.t fixed values of the 
other parameters are given in Fig. 01 We can see that all the off-forward distributions have a very strong 
i-dependence, which is similar to the i-dependence of the nuclear form factor -Fjv/aW- The nucleonic 
off- forward distribution (the leftmost panel of Fig. El nas a pronounced maximum at the point iial, 



IV. PREDICTIONS FOR PHYSICAL OBSERVABLES 



To make predictions for physical observables, we should use a model describing the internal struc- 
ture of the constituents, which is parametrized by H q u{x,£,,t) in Eq. (fTTJ) . While the resulting GPDs 
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H q l A (x, £, t) as well as all the quantities including them (cross sections, asymmetries etc.) are very 
sensitive to the details of the model used to parameterize H q u{x, £, t), we expect that the ratios of the 
nuclear-to-nucleon quantities should be less model dependent. Since modelling of the nucleon GPDs 
(especially in external meson fields) is beyond the scope of the current work, we use the simplest known 
parameterization - the double distributions 5] supplemented with the D-term |38j 

H N (x, e, t) = e l H i/n = F N(t) j i dp J J*' daS(x -/3- a£)h((3, a)q N {(3) + 9^1-^jD N (j,tj 

where ?at(x) = J2 q e q Qq/N( x ) is a standard flavour combination measured in DIS; h((3,a) — 

| (1 ^%\yf @; D N (z,t) = -E 9 e 9/ Ar / 4 ( 1 -^ 2 ) C 'i /2 (^) HI (we used Nf = 2 and neglected the D u -D d 
difference, which is formally suppressed by the 1/N C factor). For the parton distributions q q /Ki{x) we used 
the CTEQ5L parameterization |43[ . 

Since to the best of our knowledge, the fit for quark distributions in mesons was done only for 
pions ^3 but not for mesons we need, we used a simple model for the mesonic GPDs 

H v (x, Z,t)=Hf(x,Z,t) = H mes {x,U) = Y e 1 H il™° ( x > Z> *) = Js ( S ( X ~ I) ~ 6 ( X + I) ) ' (24) 

This model corresponds to a GPD of the weakly bound quark-antiquark pair. We used Nf = 2 and the 
same parameterization for both flavours since the considered mesons have isospin zero. The parameter- 
ization l|24|) has correct symmetry properties, which can be derived from the C-parity, and satisfies the 
momentum sum rule 



J y^ y xdxH q / mes (x,0,0) = 1. 



(25) 
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FIG. 4: Typical behaviour of the function R(x,£,t) (off-forward EMC effect) for 40 Ca and different values of t. 
The left panel corresponds to the parameterization 124[ , the right panel corresponds to the parameterization 1261 . 
Experimental data are from SLAC | 4q| and NMC |46j| for the forward case. 



We have neglected possible ^-dependence for the meson GPDs H mes (x,£,t) since in the kinematics of 
nuclear DVCS £ -C 1, t(r^ les ) -C 1, and give only small corrections compared to the forward case. 
Notice that this is not true for the nuclear form factors since the radii of the considered nuclei are 
(r 2 } 1 / 2 r*j 3 5 fm and, thus, values t ~ 0.1 GeV 2 are not small and the ^-dependence cannot be 
neglected. For the sake of comparison we have also investigated how our predictions change if for Hy 
and H^, we use a different model, 

H q / mes (x,^,t) = q w (x) - q w (-x) , (26) 

where q v are the pion PDFs parameterized by Owens |44| . 
The first quantity that we consider is the ratio 

which in the forward limit reduces to the ratio of the structure functions (to the leading order in as) 

R{X) = AF^ix) < (28) 

where x — Ax a- The dependence of the function R(x, £,t) on x at fixed XBj = 0.09 (HERMES) and 
different values of t is given in Fig.0] For comparison we also present the SLAC and NMC ^(| data 
for the ratio R(x). Notice that the meson contribution enhances R(x,t;,t) at small x above the pure 
nucleonic contribution previously obtained in |32t |33| . It is important to emphasize that the convolution 
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approximation, which implies the interaction with the nuclear constituents one at a time, ignores the 
physics of the simultaneous coherent interaction with several constituents, which becomes important for 
x < 0.1, see e.g. ^t}- This explains the fact that our calculation gives R(x) > 1 for x < 0.1, in dramatic 
contrast to the data. While our calculations provide a fairly good description of R(x) for x > 0.1, our 
model predicts a significant enhancement of the nuclear anti-quark distributions, which contradicts the 
experiment |48j . This problem is typical for any model of nuclear structure, which involves a significant 
fraction of non-nucleonic meson degrees of freedom |49|. 

The difference between H^ A (x, £, t) and H^ N (x, f, i) can be also probed through the DVCS beam- 
charge and beam-spin asymmetries 

where <^> is the angle between the lepton and hadron scattering planes; and a*, V denote cross-sections 
of unpolarized electron/positron and longitudinally polarized leptons, respectively [lOj. The asymmetries 
Ac(4>), Alu{4>) directly measure interference between Bethe-Heitler and DVCS amplitudes. We perform 
the evaluation of the beam-charge and beam-spin asymmetries in the kinematics similar to that used at 
HERMES Hg: {x B] ) pe rnucieon = 0.09, (Q 2 ) = 2.2GeV 2 ,(t) = -0.01 GeV 2 . In this kinematics due to 
the small value of the prefactor — t/Am 2 A , one can safely ignore the contributions of the GPD E(x,£,,t) 
and magnetic form factor i*2(£) compared to those of H(x,£,t) and F\(t), see Ref. [l(j for more details 
and a complete set of formulas describing DVCS process both for unpolarized and polarized cases. 

In HERMES kinematics, l/Q is a modest expansion parameter for the nuclear hard scattering and 
one should check if the expansion really works. For instance, due to mesons the nuclear DVCS amplitude 
squared is NOT negligible compared to the Bethe-Heitler amplitude squared (in the denominator of the 
expression for the DVCS asymmetries), though numerically it is suppressed by the factor t/Q 2 . 

Using our results for nuclear off-forward distributions, we calculate the dominant (leading-twist) 
harmonics of the beam-charge and beam-spin harmonics, 

A c c os = - f d<t> cos 4>A C {4>)- A S £B = - I dcf> sin0 A LU (<j>) . (30) 

T JO T Jo 

In order to study the role of the mesons, we give the answer for the full calculation and for the calculation, 
where the contribution of the and V mesons was neglected. The results are summarized in Table [H] 
in terms of the ratios of the nuclear to the free proton asymmetries. The second and third columns 
correspond to the calculation without the nuclear mesons; the fourth and fifth columns correspond to 
the full calculation with the meson GPDs parameterized by Eq. I)24|l : the sixth and seventh columns 
correspond to the full calculation with the meson GPDs parameterized by Eq. H2(jfl . 

One can see from Table [H] that, in the absence of the meson contributions, both asymmetries 
are practically independent of the atomic number. In this case the DVCS amplitude squared is small 
compared to the BH amplitude squared for all nuclei. 

In the presence of mesons (the full calculation) we can see that beam-charge asymmetry is a 
growing function of the atomic number whereas the beam-spin asymmetry is weakly-dependent of the 
atomic number; the DVCS amplitude increases. A least-square fit gives the following approximate A- 
dependence: Af?" A / 'A% S N cx A - 5 ; A$fi a I a lu n « A ~° ° 3 for a11 A \ thc ratio of thc DVCS amplitudes 
squared \Advcs a/Advcs n\ 2 oc A 4 29 . The natural explanation of such a behaviour is that the asymme- 
tries are large when both DVCS and BH amplitudes have comparable magnitudes and decrease when one 
of the amplitudes essentially overcomes the other. Thus in the considered kinematics the asymmetries as 
functions of the atomic number A should be small for small A (when the Bethe-Heitler amplitude domi- 
nates); they should raise when DVCS amplitude becomes comparable with the Bethe-Heitler amplitude. 
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As a consequence, the asymmetries should have their maxima. However the position of the maxima is 
strongly model dependent: the model (|24[1 predicts the maximum at moderate values of A oc 40 — 50 for 
the beam-charge asymmetry, whereas the model (|26H predicts it for A > 208 (i.e. that for all nuclei the 
A-dependence of the beam-charge asymmetry is homogeneous). 

In Fig. [5] we present the ^-dependence of the ratio of the nuclear to the free proton beam-charge 
and beam-spin asymmetries at fixed £ 0.045 (per nucleon). From Fig.JSJone can see that the qualitative 
i-dependence of both asymmetries does also change in the presence of mesons. 

It is interesting to extrapolate our results to the point t = to eliminate the "pure kinematical" 
^-dependence coming from the difference of the (r 2 ) in the nuclear formfactor -Fk(t) = Ze~( r ^*l/ 6 . The 
fit to data gives \Advcs a/Advcs n\ 2 oc A 4 - 57±0 - 17 , which is remarkably close to the value d A (0) cx A 452 
obtained in Section^] In the asymmetries nuclear formfactors contract and we get for t — approximately 
the same power as for finite t. 



Nucleus 


A COS / A COS 

A C Al A C N 


A sin / A s in 
A LU AI A LU N 


A COS 1 A COS 

A C A/ A C n 


a sin j a sin 
A LU A/ A LU N 


A COS 1 A COS 

A C AI A C N 


A sin / A sin 
A LU Al A LU N 


12 C 


2.45 


1.85 


4.61 


2.49 


1.573 


2.222 


16 Q 


2.43 


1.83 


5.41 


2.33 


1.905 


2.270 


40 Ca 


2.38 


1.79 


7.34 


1.60 


3.276 


2.180 


90 Zr 


2.59 


1.93 


6.80 


0.81 


4.879 


2.104 


208 Pb 


2.42 


1.07 


6.12 


0.31 


6.077 


0.998 



TABLE II: The ratios of the nuclear to the free proton asymmetries for different nuclei. The second and third 
columns correspond to the calculation without the nuclear mesons; the fourth and fifth columns correspond to 
the full calculation including the meson contribution with the internal structure model given by (1241 : the sixth 
and the seventh columns correspond to the parameterization 12611 . 



It is interesting to compare our predictions for the DVCS asymmetries to the HERMES preliminary 
result on Neon 19]. Taking the ratio of the nuclear to the proton single-spin asymmetries measured by 
HERMES, one readily finds 

Asin 

LU A = 1.22 ±0.26, (31) 

A sin ^ ' 

A LU N 

where we have added the proton and neon experimental errors in quadrature. In the similar kinematics, 
the linear interpolation of our results from Table ITU gives 

-ffA«2.1±0.2. (32) 

A LU N 

A possible explanation of the discrepancy between the experimental result and the theoretical prediction 
(a simple combinatoric counting also suggests the enhancement) was suggested in p7) . The explanation 
consists in the observation that, depending on t and the experimental resolution, the contribution of 
incoherent nuclear scattering might significantly reduce the nuclear DVCS asymmetries. 



V. CONCLUSIONS 



We performed a microscopic evaluation of nuclear GPDs for spin-0 nuclei in the framework of the 
Walecka model. We found that the meson (non-nucleonic) degrees of freedom of the Walecka model play 
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0.05 0.1 -I 

0.05 0.1 -1 

FIG. 5: The ratio of the nuclear to the free proton beam-charge and beam-spin asymmetries as functions of the 
momentum transfer squared t at fixed £ « 0.045 (per nucleon) for the C 12 nucleus. The results are presented 
for the full calculation with the parameterization 12(il (solid curves) and for the calculation ignoring the meson 
contribution (dashed curves). 



a dramatic role in nuclear GPDs. 

We observed a non-trivial A-dependence of the first moment of the nuclear D-term, c?a(0) oc A 2 26 , 
which confirmed the prediction of M. Polyakov made in the framework of the nuclear liquid drop 
model 20]. This result demonstrated that contrary to the assumptions of the random phase approx- 
imation, mesonic degrees of freedom dominate (Ia{0) and cannot be neglected. 

Using the resulting nuclear GPDs, we studied the A and ^-dependence of the beam-charge and 
beam-spin DVCS asymmetries in the HERMES kinematics. We found that due to the mesonic degrees of 
freedom, the nuclear DVCS amplitude squared has a rapid A-dependence, |A£>ycs| 2 oc A 4 - 29 , and becomes 
much larger than the Bethe-Heitler amplitude squared. This dramatically affects the A-dependence of 
the asymmetries. We found that A C ° S A /A C ° S N oc A - 5 ; A s ™ A /A s ™ N oc A" 03 for all A. This behaviour 
should be compared to the case, when the meson degrees of freedom are neglected. In this case, the 
DVCS asymmetries are virtually A- independent. 
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APPENDIX: DERIVATION OF THE GPDS 



In this section we establish the formulas used for evaluation of GPDs in the functional integrals 
approach |50| . The generalized parton distributions are defined as 

H +( x >Z> t )=2^P+ J ^ XP+Z ~{P'\d+cl>(-~n-)d + cj>(^)\P), 



^ e **P+*-W V+ (-in-)V+(-n-)\P). (A.l) 
2xP+ J 2ir 2 2 y ' 

For the sake of brevity we denote <&(x) = {<f>(x), V(x), &(x), ^(x)} and notice that in all three cases we 
have to evaluate the general matrix element 

(P>\* (1\ * (_ 1\ IPS = I m*)DHx)Dn*)DV{x)* (f) * H) z %S (A 9] 

V 1 \2/ V 27 1 ' fD<t>(x)D*(x)D*(x)DV(x)e iS ' 1 ' 

which, up to inessential kincmatical factors and integration over the light-cone separation z, will give us 
the GPDs dj). 

By definition we should integrate in 1A.2|) over all configurations which satisfy the asymptotic 
conditions 50] 



lim t) » $ p < p ' (x, t) , (A.3) 

t— >±oo 



or in more explicit form 



lim w $ a (£- V'fi), lim t) w $ s (x - Vt 2 ) ■ (A A) 

ti — > + 00 t2 — > — OO 

We reparameterize the field configurations as 

-> $ a (x-Jt(t))+5$(x-j?(t)), 
L>$ -► / L><S$L>A(i) , (A.5) 



where the integral over the zero mode - position of the soliton X(t)- should be taken exactly, while 
the integral over the shifts 6<fr(x) may be evaluated in the saddle-point approximation. Notice that the 
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standard mean field approximation ignores all the loop corrections and thus implicitly uses the smallncss 
of the coupling constant g. In this case the evaluation of a simple path integral for the center of mass 
motion gives 

J D *(t)* (|) * (--) z lS « ^ J dX x J dX 2 <S> (| - X X ) $ (-| - J? 2 , x 



t>z°/2 

const 



Yl dX(t)e iM ' x ^/ 2 



dX x dX 2 exp 



f<-z°/2 

»M g (Xi - X 2 ) 2 
2^o 



-z /2<i<z /2 



exp( l (P'X 1 - ?i 2 ))$ s (| - Xi) $ s (-| - X 2 )a.6) 



where Xl = J?(f = z°/2);X 2 = X(f = — z°/2) and continuity of the path as well as proper boundary 
conditions 

*(tl) = ft, X(t 2 ) = q 2 , lim fi=y', lim ^ = F 

are implied; the const represents (divergent) kinematical factors which will contract with the same factors 
in the denominator of ljA.21) . From the definition 1A.1|) we can see that the most essential are small 



z r~j 1/P+. The term const exp 



2z 



Z^/zq is not analytical at the point z$ = and is strongly 



suppressed for zq ^ 0. Hence we may safely replace it with its limit 



lim exp 



iM s {X 1 -X 2 ) 2 
2z 



S(X 1 - X 2 ) 



Thus finally we arrive to the well-known result [FH l52l | 

(^(|)*(-f)i p '^/ to -(|-^)*- 



- 2 -X\e 



■AX 



(A.7) 



which is valid in the leading order in the coupling constant: to evaluate the matrix element, one should 
replace all the fields with the solutions of the Hartree-Fock equations centered at the point X and integrate 
over the center of mass. 
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